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Most investigations devoted to the conditions for adiabatic quantum computing are based on the 
first-order correction (*I'ground(t)|i?(t)|^'cxcitcd(t))/AiJ^(t) <C 1. However, it is demonstrated that 
this first-order correction does not yield a good estimate for the computational error. Therefore, a 
more general criterion is proposed, which includes higher-order corrections as well and shows that 
the computational error can be made exponentially small - which facilitates significantly shorter 
evolution times than the above first-order estimate in certain situations. Based on this criterion and 
rather general arguments and assumptions, it can be demonstrated that a run-time T of order of 
the inverse minimum energy gap AiJmin is sufficient and necessary, i.e., T = 0{/S.E'^^. For some 
examples, these analytical investigations are confirmed by numerical simulations. 

PACS numbers: 03.67.Lx, 03.67.-a. 



I. INTRODUCTION 

With the emergence of the first quantum algorithms, 
it turned out that quantum computers are in principle 
much better suited to solving certain classes of problems 
than classical computers. Prominent examples are Shor's 
algorithm jjyj for the factorization of large numbers into 
their prime factors in polynomial time and Grover's al- 
gorithm I3| for searching an unsorted database with N 
items reducing the computational complexity from the 
classical value 0{N) to 0{^/N) on a quantum computer. 

Unfortunately, the actual realization of usual sequen- 
tial quantum algorithms (where a sequence of quantum 
ates is applied to some initial quantum state, see, e.g., 
) goes along with the problem that errors accumu- 
late over many operations and the resulting decoherence 
tends to destroy the fragile quantum features needed for 
the computation. Therefore, an alternative scheme has 
been suggested 0, where the solution to a problem is 
encoded in the (unknown) ground state of a (known) 
Hamiltonian. By starting with an initial Hamiltonian 
with a known ground state and slowly evolving to the fi- 
nal Hamiltonian Hi with the unknown ground state, e.g., 
H{t) = [1 — s{t)\Hi + s{t)Hf, adiabatic quantum comput- 
ing makes use of the adiabatic theorem which states that 
a system will remain near its ground state if the evolu- 
tion s(t) is slow enough. Since there is evidence that the 
ground state is more robust against decoherence 0, 0, j 
this scheme offers fundamental advantages compared to 
sequential quantum algorithms. 

However, determining the achievable speed-up of adi- 
abatic quantum algorithms (compared to classical meth- 
ods) for many problems is still a matter of investigation 
and debate, see, e.g., |1 11 [iSM 111 111 111 . For ex- 
ample, it has been argued in |10| that all conventional 
(sequential) quantum algorithms can be realized as adi- 
abatic quantum computation schemes with polynomial 
overhead via the history interpolation (polynomial equiv- 
alence). For an adiabatic version of Grover's algorithm, a 
constant velocity s implies a linear scaling of the run-time 



T — 0{N), whereas a suitably adapted time-dependence 
s(t) yields the known quadratic speed-up T = 0{^/N), 
cf. [la. . Whether adiabatic algorithms of NP com- 
plete problems such as 3-SAT can be even more efficient 
than this quadratic speed-up is still not clear, see, e.g.. 

Si. 

In this paper, we derive a general error estimate as a 
function of the run-time T (the main measure for the 
computational complexity of adiabatic quantum algo- 
rithms) for very general gap structures Ai?(s) and in- 
terpolation velocities s{t). 



II. ADIABATIC EXPANSION 

The evolution of a system state |^'(t)) subject to a 
time-dependent Hamiltonian H{t) is described by the 
Schrodinger equation {h = 1) 



t\^{t))^H{m{t)). 



(1) 



Using the instantaneous energy eigenbasis defined by 
H{t)\n{t)) = En{t)\n(t)), the system state {^(t)) can be 
expanded to yield 



1^(0) = J2 «»(^) I J En{t')dt' \ \n{t)) . (2) 



Insertion into the Schrodinger equation yields - after 
some algebra - the evolution equations for the coefficients 
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dt 



{m\H\n) 



xexpi -2 j AE^rn{t')dt' [ (3) 
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with the energy gap AEnm (t) = En (t) — Em (t) and the 
Berry phase |l5l | 



7„(i)=* J dt' {n{t')\n{t')) . 





(4) 



If the external time-dependence H is slow (adiabatic evo- 
lution), the right-hand side of Eq. Q is small and the 
solution can be obtained perturbatively. After an inte- 
gration by parts, the first-order contribution yields 



an- 



AEL 



where ipnm G R denotes a pure phase. Consequently, if 
the local adiabatic condition 



{m\H\n) 
A£;2 



£ < 1 



(6) 



is fulfilled for all times, the system approximately stays 
in its instantaneous eigen (e.g., ground) state through- 
out the (adiabatic) evolution. This above constraint has 
frequently been used as a condition for adiabatic quan- 
tum computation However, since the solution to a 
problem is encoded in the ground state of the final Hamil- 
tonian in adiabatic quantum computation schemes, it is 
not really necessary to be in the instantaneous ground 
state during the dynamics - the essential point is to ob- 
tain the desired ground state after the evolution. Since 
the external time-dependence H could realistically be ex- 
tremely small (or even practically vanish) at the end of 
the computation t = T, the first-order result lO does 
not always provide a good error estimate. Similar to 
the theory of quantum fields in curved space-times [T^ . 
the difference between the adiabatic and the instanta- 
neous vacuum should not be confused with real excita- 
tions (particle creation). Therefore, it is necessary to go 
beyond the first-order result above and to estimate the 
higher-order contributions. 



III. ANALYTIC CONTINUATION 

Evidently, the Schrodinger equation is covariant under 
simultaneous transformations of time and energy, such 
that the runtime of any adiabatic algorithm can be re- 
duced to constant if the energy of the system is modified 
accordingly |13|. Here we want to exclude a mixing of 
these effects and will therefore assume 



Ti{H[s{t)]} = const. V s e [0, 1] 



(7) 



where < s(i) < 1 is an interpolation function which 
will be specified below. In practice, the above condition 
can even be relaxed to the demand that the trace should 
not vary by orders of magnitude (during < s < 1). 



With suitable initial and final Hamiltonians Hi and H{, 
the above condition can be satisfied for all s by using the 
linear interpolation scheme 



H{t) = [1 - s{t)] Hi + s{t)Hi 



(8) 



but other schemes are also possible (see sectional. For 
simplicity, we restrict our considerations in this section 
to a non-degenerate (instantaneous) ground state n = 
and one single first exited state m = 1 with AE = AEio . 
(Multiple excited states will be discussed in section 0) 
Similarly, all energies will be normalized in units of a 
typical energy scale corresponding to the initial/final gap, 
A£;(0) = 0(1) and A£:(l) = 0(1). We classify the 



(5) dynamics of s(t) via a function h{s) > 



^^^AEis)his), 



(9) 



where the function h{s) > is constrained by the condi- 
tions s(0) = and s{T) = 1. Insertion of this ansatz into 
Eq. Q yields the exact formal expression for the non- 
adiabatic corrections to a system starting in the ground 
state, i. e., with ai(0) = one obtains after time T 



ai(l)e 



-i7i(l) 



ds ao(s)e 



Foijs) 
AEis) 



X exp ■ 



ds' 
W) \ ' 



(10) 



with the matrix elements Fnm{s) — ('Ti(s)|-ff'(s)|n(s)) 
which simplify in the case (jS]) of linear interpolation to 
Fnm{s) — {m{s)\{H{ — Hi)\n{s)). The advantage of the 
form in Eqs. © and IjlOf) lies in the fact that different 
time-dependences s{t) and hence different choices for h{s) 
solely modify the exponent. 

We assume that all involved functions can be analyt- 
ically continued into the complex s-plane and are well- 
behaved near the real s-axis. Given this assumption, we 
may estimate the integral in Eq. HlOfl via deforming the 
integration contour into the lower complex half-plane (to 
obtain a negative exponent - which is the usual procedure 
in such estimates) until we hit a saddle point, a singu- 
larity, or a branch cut, see Fig. ^ Deforming the inte- 
gration contour into the upper complex half-plane would 
of course not change the result, but there the integrand 
is exponentially large and strongly oscillating such that 
the integral is hard to estimate. Since the gap AE{s) 
usually has a pronounced minimum at Smin G (0, 1), the 
first obstacle we encounter 17] will be a singularity at s 
close to the real axis, i.e., |3(s)| <C 1 and 3fJ(s) w Smin, 
where AE{s) = 0. 

Let us first consider a constant function h{s) ~ h: 
Assuming h ^ I (i.e., slow evolution), the exponent in 
Eq. ((Tn|l acquires a large negative real part for 3(s) < 
and thus the absolute value of the integrand decays 
rapidly if we depart from the real s-axis in the lower com- 
plex half-plane. Imposing the even stronger constraint 
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FIG. 1: [Color Online] The original integration contour (black 
line along real axis) of equation IIUII is shifted to the complex 
plane (curved line). The gap structure AE{s) leads to singu- 
larities near the real axis [green hollow circles, here displayed 
for 2a = 4 in Eq. Il7ll 1 , which limit the deformation of the in- 
tegration contour. The integral in the exponent (dashed line) 
in equation IIUII ranges from to s', which gives rise to a real 
contribution to the exponent off the real axis only. 



h <C 13^(5)1 <C 1, the decay of the exponent dominates 
all the other s-dependences [71 (s), Foi(s), and AE{s)] 
since their typical (minimum |l7| ) scale of variation is 
|S(s)| ^1. In view of the complex continuation of 
Eq. the same applies to the amplitude ao(s). As a 
result, the above integral ()10|l will be exponentially sup- 
pressed ~ exp{-0(|5(s)| /h)} if /i< |3(s)| < 1 holds, 
which (as one would expect) implies a large evolution 
time T via the side condition s(T) = 1. 

The general situation with varying h(s) can be treated 
in complete analogy - the integral in Eq. (|10|l is sup- 
pressed provided that the condition 



h{0) + h{l) < 1 A 3? 



/ SR(5)-l-ia(s)/2 



V 




»1 (11) 



holds for all singularities s (and saddle points etc.) in the 
lower complex half-plane (which determine the deforma- 
tion of the integration contour) . Together with 



T 



ds 



AE{s)h{s) 



(12) 



this determines an upper bound for the necessary runtime 
T of the quantum adiabatic algorithm. 

Note that the constraint s ^ AE derived from 

h <C |3^(s)| is not necessarily equivalent to s ^ AE'^, 
which one would naively deduce from Eq. (jSJ. 



IV. EVOLUTION TIME 

The general criterion in Eq. I|ll|) can now be used to 
estimate the necessary run-time via Eq. (fT^ . Typically, 
the inverse energy gap l/AE{s) is strongly peaked (along 
the real axis) around 5i(s) with a width ^3 of order 
|Q'(s)|. Therefore, assuming h{s) to be roughly constant 
across the peak and respecting h |pcak^ yields 



the following estimate of the integral in Eq. (|12|l 



T = 0{AE-2,) , 



(13) 



where Ai^min denotes the minimum energy gap. Note 
that this estimate is only valid for one (or a few) relevant 
excited state (s) - multiple excited states will be discussed 
in section Ivl 

Intuitively, the same order of magnitude estimate for 
the evolution time can also be derived from the local 
adiabatic condition ® : Inverting this condition, wc find 
the relationship 



T 



1 



ds 



Foijs) 
A£;2(s) 



(14) 



Assuming that i^oi(s) does not oscillate strongly, e.g., 
that the ground state of H{s) travels on a reasonably 
direct path from the initial to the final state, we can 
make the following estimate 



T = 



OiAE-l) 



ds 



Fpiis) 
AE{s) 



(15) 



Now we may exploit the advantage of the representa- 
tion in Eq. H1U|I . which is valid for general dynamics s{t) 
corresponding to different functions h{s) and hence for 
arbitrary evolution times T. In the limit of very fast 
evolution T —* (which implies h 00), we have large 
excitations ai{T) = 0{1) and thus the remaining inte- 
gral in the above equation can be estimated via inserting 
this limit into Eq. l(TU|l : 



ds 



f 01(g) 
AE{s) 



Oil). 



(16) 



By comparing Eqs. 1)16(1 and H14I) . we again obtain the es- 
timate (|13|l . Note that the quantities Foi(s) and AE{s) 
appearing in the integrals in Eqs. H14I16|I do not depend 
on the dynamics s{t) which allows us to perform the in- 
tegration independently of s{t). 



A. Gap Structure 

Let us illustrate the above considerations by means of 
the rather general ansatz for the behavior of the gap 



AE{s)^[is- 



1/6 



(17) 



with the minimal gap < AE'min ^ 1 at Smin G (0, 1), 
6 > 0, and a E N+. An avoided level crossing in 
an effectively two-dimensional subspace corresponds to 
2a = b — 2. This is the typical situation if the commu- 
tator of the initial and the final Hamiltonian [Hi, Hf] is 
small, since, in this case, the two operators can almost 
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FIG. 2: [Color Online] Runtime scaling of the adiabatic 
Grover search for different interpolation functions s{t) and a 
target fidelity of 3/4. Solid lines represent fits to full symbol 
data for A'^ > 100 and shaded regions correspond to fit un- 
certainties (99% confidence level). These uncertainties arise 
from the finite resolution when determining the necessary run- 
time. Hollow circles represent calculations with smoothed 
C°°-interpolations (compare dotted lines in figure 2] and sec- 
tion EJi whereas hollow boxes correspond to the nonlinear 
interpolation example in section IVI 



be diagonalized independently and thus the energy levels 
are are nearly straight lines except at the avoided level 
crossing(s), where [Hi, H{] becomes important. In the 
continuum limit, such an (Landau-Zener type) avoided 
level crossing corresponds to a second-order quantum 
phase transition. The finite-size analogue of a third- 
order phase transition corresponds to a = & (and ac- 
cordingly for even higher orders), which may occur if 
[Hi, H{] is not small or if the interpolation is not linear, 
i.e., H{s) ^ [l-s]Hi + sHf. 

The inverse gap 1/ AE{s) has singularities around Smin 
at 3(s) = 0{AE'^^^), compare Fig.^ The total running 
time T for different choices of h{s) = adAE'^{s) satisfying 
the criterion Ijlll) can be obtained from Eq. H12|l . Here, 
the exponent d determines the scaling of the interpolation 
dynamics, whereas the coefficient ad is adapted such that 
s(T) = 1, cf. Eqs. © and (P|l . 

For 2a{d+l)/b>l one easily shows that 
l/ad = OiTAE'^l'^^^") satisfies the criterion HH) 
with the evolution time obeying T — 0{AE^^l^J. If d is 
smaller, the necessary evolution time will be larger. In 
Table HJ the scaling of the run-time (for two examples of 
the gap structure) is derived for three cases: 

a) constant velocity s — a_i, i.e., d = —1, 

b) constant function h{s) = ap, i.e., d — 0, and 

c) the local adiabatic evolution with h{s) — aiAE{s), 
i.e., d = +1, investigated in p^ . 



AE{s) = 



d=-l 
d = 
d>l 



,/is-l/2y + AE-i,^ 



AEZl 



^is-l/2r + AE- 



AE~^/^ 

mm 
-1 
mill 



AE, 
AE- 



TABLE I: Scaling of the runtime T necessary to obtain a 
fixed fidelity for different gap structures (top row) and varying 
interpolation velocities (first column). The best improvement 
possible scales as the inverse of the minimum gap AE^^. 



B. Grover 's Algorithm 

In the frequently studied adiabatic realization of 
Grover's algorithm (see, e.g., [l^ ITsL IT^ ) the initial 
Hamiltonian reads Hi — \— |in)(in| with the initial su- 
perposition state \va) — ^^^^ \x) / \fN , and the final 
Hamiltonian is given by iJf = 1 — |?i')(w|, where |w) de- 
notes the marked state. In this case, the commutator 
is very small [Hi, Hi] — (|in)(w| — \w) {m\)/\/N and one 
obtains for the time-dependent gap |l4| 



AE{s) 




(18) 



Comparing with Eq. H17|l , we identify Ai^inin ~ 1/ yN 
and 2a = 6 = 2 (the pre-factor does not affect the scaling 
behavior). Consequently, our analytical estimate implies 
T = 0{N) for d = -1, T = 0{^/N\nAN) for d = 0, and 
T = oI^/n) for d > 0. 

We have solved the Schrodinger equation numerically 
by using a fourth order Runge-Kutta integration scheme 
with an adaptive step-size [ig. By restarting the code 
with different T until agreement with desired fidelity was 
sufficient, we could confirm these runtime scaling pre- 
dictions numerically, see Fig. [2| The dependence of the 
final error on the run-time T for fixed N — 100 and con- 
stant h is depicted in Fig. O where the exponential de- 
cay becomes evident. The evolution of the instantaneous 
ground-state occupation is plotted in Fig. 0] for the three 
different dynamics. 



V. FURTHER GENERALIZATIONS 
A. Adiabatic Switching 

From an experimental point of view, the time- 
dependence of the Hamiltonian will most certainly vanish 
asymptotically H{t < 0) = H{t > T) = or at least be 
neghgible - which automatically imphes /i(0) = h{l) = 0. 
Furthermore, realistic Hamiltonians should be described 
by C°°-interpolations [Natura non facit saltus). 

By using a C°°-test function which was matched at 
ti = O.IT and ^2 = O.QT to the usual dynamics s(t) 
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FIG. 3; [Color Online] Final error probability |ai(r)^| as 
a function of run-time T for Grover's algorithm with A'^ = 
100 and h — const. The oscillations stem from the time- 
dependence of ao in Eq. . The solid (blue) line represents 
the second-order perturbative solution of Eq. IjlO^ . 



(compare dotted lines in figure0]bottom panel), we have 
implemented an interpolation scheme with such an adia- 
batic switching on and off s(0) — s{T) = 0. For the in- 
vestigated adiabatic implementation of the Grover search 
routine, this scheme does not affect the final result con- 
siderably. The reason for this robustness lies in the fact 
that the matrix element Fnm is peaked around s — 1/2 
and h{0) as well as h{l) are small enough already without 
the adiabatic switching on and off. Therefore, one can 
expect the dominant non-adiabatic corrections to arise 
from the behavior around the minimum gap, which was 
unaffected by the test function. This is also confirmed by 
the scaling of the runtime versus the system size, compare 
the hollow circle symbols in figure |21 which is basically 
unchanged. 

However, the situation is completely different for the 
example considered in section FV CI below. There, the ex- 
ponential suppression of the final error as a function of 
the run-time requires a smooth C°°-interpolation - with 
other dynamics such as C" (just continuous) or (dif- 
ferentiable once), the final error is merely polynomially 
small, cf. figure |31 



B. Nonlinear Interpolation 

Although we have chosen a linear interpolation 
scheme (jSJ in order to satisfy the trace constraint ([TJl , the 
presented analysis can be generalized easily to more gen- 
eral non-linear interpolations. [Note that, linear refers 
to the straight connection line between initial and final 
Hamiltonian in equation Q and should not be confused 
with the different velocities s{t) at which this line is tra- 
versed.] The argumentation based on the analytic con- 
tinuation works in the same way provided that the func- 
tional dependence Hni{s) = f{Hi,H{,s) does not involve 
extremely large or small numbers. 




FIG. 4: [Color Online] Evolution of the interpolation function 
s{t) (bottom panel), the spectrum a[s{t)] (middle panel), 
and the occupation of the instantaneous ground state (top 
panel) versus the rescaled time r — t/T for an adiabatic 
Grover search problem with A'' = 100 states. For each inter- 
polation (different line styles) , T was adapted to reach 99% of 
final fidelity. Thin dotted lines represent C°°-interpolations 
smoothed with a test function. 



As an illustrative example, we consider the Grover 
search with the same initial and final Hamiltonians but 
a quadratic interpolation scheme 



i?„i(s) = [(l-s)Hi + sHf]^ + s{l-s) 



2N 



iV2 



(1 - sfHi + s^Hi 



+s(l-s) 



{Hi,Ht} 



2N-2 



(19) 
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FIG. 5: [Color Online] Evolution of the final and the maxi- 
mum intermediate (red line) excitations with the runtime T 
for the example l|2()^ . The exponential falloff in the final ex- 
citations is only visible, if a smooth C°° -interpolation (black 
circles) is used, whereas the scaling of the intermediary ex- 
citations (red line) is always polynomial. The suppression of 
the final error for C° or C'^-interpolations (blue squares and 
green crosses) is also merely polynomial. 



where {•,•} denotes the anti-commutator. The identity 
operator 1 has been added in order to ensure Tr{_ff,ii} — 
N ~ 1, cf. equation Q. Although the spectrum of this 
non-hnear interpolation is slightly distorted compared to 
the linear one, the fundamental gap is the same as in 
equation (|18|l . and hence same interpolation functions 
s{t), applied to the above Hamiltonian, should reproduce 
the aforementioned scaling predictions. This is confirmed 
by the numerical analysis of the scaling behavior - the 
results of the non-linear interpolation are basically indis- 
tinguishable from those of the previous example (linear 
interpolation), compare the hollow box symbols in fig- 
ure |2 



C. Degeneracy 

So far, we have restricted our considerations to the in- 
stantaneous ground state and a single first excited state. 
Let us now consider a very simple example (see also 0) 
in which there is still a unique ground state, but many 
degenerate first excited states: In terms of single-qubit 
Pauli matrices ax and CTz , the M-qubit Hamiltonian reads 



H{s) 



M 



i(i) 



(20) 



where we have used a linear interpolation © for sim- 
plicity. In this example, the Hamiltonian can be decom- 
posed completely into independent and equal single-qubit 
contributions and hence the time-evolution operator fac- 
torizes, i.e., it is sufficient to solve the dynamics of a 
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scale factor T=t/T 

FIG. 6: [Color Online] Occupation of the instantaneous 
ground state and some selected computational basis states 
for the Hamiltonian in (1201 for an M = 8 qubit system. Tem- 
porarily, the system leaves the instantaneous ground state, 
but the runtime T has been adjusted such that the final fi- 
delity is 99%. 



single qubit. Furthermore, the Hamiltonian is invari- 
ant under any permutation of the qubits. The instan- 
taneous ground states for all values of s are symmetric 
under this permutation group and hence unique, but the 
first excited states are not - leading to a Af-fold degen- 
eracy (i.e., there are M equivalent first excited states). 
Hence, the fundamental gap between the ground state 
and each one of these first excited states is the same as 
for one qubit and thus independent of the number of 
qubits AE{s) = y^l - 2s(l - s). 

In some sense, this simple example represents a lim- 
iting case opposite to Grover's algorithm: The energy 
gap AE{s) and the matrix elements Fnm{s) do not scale 
with the number M of qubits and the Fnm are neither 
small initially nor finally. Instead, the scaling with sys- 
tem size manifests itself in the Af-fold degeneracy of the 
first excited states. As a result of the M-independent gap 
structure, the adiabatic switching is crucial for achieving 
the exponential suppression of the final error. Figure [51 
displays the final error probabilities for a smooth C°°- 
interpolation and for and C^-interpolations for com- 
parison. These numerical simulations confirm that the 
falloff is exponential in the C°°-case but merely polyno- 
mial for C° and 

Another interesting point of this simple example is the 
difference between the intermediate and the final occupa- 
tion of the ground state, see figures [HI and [5| According 
to the first-order result in Eq. ||3J) and the aforementioned 
factorization of the time-evolution operator, the interme- 
diate excitation probability scales as 



Pint 



El 

m>0 



M 



(21) 



since the gap /S.E is independent of M. On the 
other hand, the final error probability (assuming a C°°- 
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interpolation) is exponentially suppressed 

Pfin^O{Mexj>{-TAE}) ^0{Mexp{-T}) , (22) 

and hence the two error probabilities can be vastly dif- 
ferent Pint ^ Pfinj cf. figure IHI In fact, by increasing the 
number of qubits, the occupancy of the instantaneous 
ground state can be made arbitrarily small. Moreover, 
the run-time condition derived from the first-order result 
in Eqs. © and (jHl 

To = 0{Vm) , (23) 

yields a scaling which is far too pessimistic compared 
with the correct final error probability assuming a C°°- 
interpolation 

Too-O(lnM). (24) 

Note that non-smooth interpolations (e.g., C*' or C^) 
would also yield a polynomial scaling T — 0{M'-^) sim- 
ilar to Eq. (|23|l . On the other hand, the scaling behav- 
ior in Eqs. H22I) and 1)24(1 is just what one would ob- 
tain by immersing the system in Eq. (|20|l into a zero- 
temperature environment and letting it decay towards 
its ground state. Therefore, using non-smooth interpo- 
lations (e.g., C° or C^) or naively demanding the first- 
order estimate in Eq. ©, the adiabatic algorithm would 
be even slower than this simple decay mechanism. 

VI. SUMMARY 

The instantaneous occupation of the first excited state 
during the adiabatic evolution in Eqs. jSJ) and © does 
not provide a good error estimate. Instead, a better esti- 
mate is given by the remaining real excitations after the 
dynamics. For the example plotted in Fig.^ the instan- 
taneous excitation probability exceeds 10% at intermedi- 
ate times - whereas the final value is 1%. This is even 
more drastic for the example in section FV CI see figure |S| 
where the two values and hence the inferred run-times 
can differ by orders of magnitude. 

Moreover, the final error can be made extremely - in 
fact, with h{0) + h{l) <^ 1, exponentially - small 

ai(r) = O (^h{0) + h{l) + cxp , (25) 

cf. Fig. 13 For the Grover example, the last term was 
dominant, whereas in the general case the smallness of 
the first two terms can be ensured by using smoothed 
C°°-interpolations, i.e., adiabatic switching - which is a 
more realistic ansatz anyway. 

Based on general arguments, the optimal run-time (in 
the absence of degeneracy, cf. section IV C|l scales as 
T = 0{AE^l^) contrary to what one might expect from 
the Landau-Zener formula (with T oc AE^f^). In 
view of the fact that the minimum energy gap Ai?min is 



a measure of the coupling between the known initial state 
and the unknown final state, this result is very natural. 

For the Grover algorithm, it is known that the 
ViV-scaling is optimal This optimal scaling 

T — 0{AE^l^) can already be achieved with interpola- 
tion functions s{t) which vary less strongly (e.g., d — 0) 
than demanded by locally fl^ adiabatic evolution {d — 1) 
- and hence should be easier to realize experimentally. 

Unfortunately, a constant velocity with d = —1 does 
not produce the optimal result in general. The Grover 
example has the advantage that the spectrum can be de- 
termined analytically, which is for example not the case 
for the more involved satisfiability problems There- 
fore, some knowledge of the spectral properties AE{s) is 
necessary for achieving the optimal result T = 0{AE^^l^^) 
also in the general case of adiabatic quantum computing. 
For systems with an analytically unknown gap structure, 
some knowledge about the spectrum can be obtained by 
extrapolating the scaling behavior of small systems. 

A related interesting point is the impact of the gap 
structure (corresponding to 2"*^ or 3"^ order transition 
etc.) in Eq. (|17|) . The derived constraint for the velocity 
at the transition s ^ |5(s)| AE is only for 2"'^-order tran- 
sitions equivalent to s ^ AE'^, which one would naively 
deduce from Eq. ©. 

Note that the improvement T — 0{AE^l^^) compared 
with the conventional runtime estimate T = 0{AE^f^) 
is merely polynomial (same complexity class). Though 
this is not as impressive as an exponential speedup, in 
practice a polynomial improvement may be useful. For 
time-dependent Hamiltonians where the inverse of the 
minimum gap scales exponentially with the size of the 
problem, we would still expect an exponential scaling of 
the runtime T required to reach a fixed fidelity (as in 
section llV Bp . On the other hand, the exponential sup- 
pression of the final error in Eq. H25|) may become impor- 
tant in certain cases such as in the presence of degeneracy 
and may well yield an exponential speedup in comparison 
with the conventional estimate, see section IV CI 

In some sense, the two examples in sections II V Bl 
and IV Cl represent two simple extremal examples for adi- 
abatic quantum computing regarding the scaling of the 
gap and the degeneracy. For more complicated situations 
such as satisfiability problems "ij, both properties have 
to be taken into account simultaneously. 
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Note added than-polynomial decrease of the final error ai(T), have 

been demonstrated rigorously for a class of Hamiltonians 
Recently, two of the main results of this article, i.e. , the using methods of spectral analysis 13 . 
optimal run-time scaling T = 0{/!^E^_^^ and the faster- 
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